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Generating multivariate extreme value distributions 



Helena Ferreira Department of Mathematics, University of Beira Interior, Covilha, Portugal 

Abstract: We define in a probabilistic way a parametric family of multivariate extreme value distribu- 
tions. We derive its copula, which is a mixture of several complete dependent copulas and total independent 
copulas, and the bivariate tail dependence and extremal coefficients. Based on the obtained results for 
these coefficients, we propose a method to built multivariate extreme value distributions with prescribed 
tail/extremal coefficients. We illustrate the results with examples of simulation of these distributions. 
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CN 1 Introduction 

u 

The construction of multivariate distributions has its motivation in Probability Theory, Biostatis- 
tics, Economics, and there is no need to demonstrate its importance nowadays. In particular, 
multivariate extreme value (MEV) distributions provide models for joint extreme events which 
require special care from practitioners since rare events can have serious natural and economic 

i— i impact. 

The MEV distributions composes an important class of positive dependent distributions with de- 

0_| pendence in the extreme upper tails, with this property known as tail dependence ( J9lfl4]). A 

misleading evaluation of the dependence in the tails of a multivariate distribution may lead to an 
underestimation of risks from crises events ( [2]). 

The MEV distributions can be constructed through methods such as mixtures, stochastic repre- 
sentations and limits ( J9l). They also can be obtained via copulas and suitable techniques of 
extra-parametrisation ( [6J[7 11 12]). 



y—i The dependence structure of a MEV distribution is completely characterised by its dependence 

function ( |14| ). However this function cannot be easily inferred from data and simple dependence 
measures are welcome. 

The most popular of the dependence measures for bivariate distributions are the upper tail de- 
pendence coefficient A and the extremal coefficient e, introduced far back in the sixties ( fl8||l9] ), 

which determines each other in the extreme value distributions by the relation e = 2 — A. 
C^) ... 

Our paper is a contribution for the main issue of the development of higher dimensional models 

capturing dependencies. This issue covers quite a large spectrum of techniques and applications 
y—i where copulas become ladies of quality since they are able to yield any kind of dependence struc- 

ture independently of marginal distributions. An application of copula functions is the simulation 
of multivariate distributions with dependent observations. By departing of the right copula or 
a mixture of copulas any dependence structure may be reached. For a review of different esti- 
mators, some model selection tests for copula functions and methods of its simulation see, for 
instance, |4[|8| |13| . However, in general copulas cannot be characterised in a simple way that pro- 
vides a simple algorithm to simulate data. 

How to get out of the way of selection of a candidate copula, its estimation and simulation of 
dependent random observations? One proposal is to estimate non-parametricaly the bivariate ex- 
tremal dependences ( |5| |16| ) and to construct in a probabilistic way a multivariate extreme value 
distribution with these prescribed extremal dependences, which can be used to generate random 
observations. 

Here we propose a method to built a MEV distribution with prescribed tail/extremal dependence 
coefficients. 

To the best of our knowledge, the problem was treated by [I] for elliptical distributions and by [3] 
for the case where each bivariate tail dependence coefficient A,j between the margins Xi and Xj 



has the representation A^ = 2 — (1 + a; — a/) 1 '", where < a\ < ... < ad < 1 < a and d is the 
dimension of the random vector. 

Our solution is of stochastic representation kind and can be used in practical problems where one 
needs to buil a stochastic model in a situation where the extremal dependence is known or esti- 
mated and the knowledge or estimates of marginal distributions are available. 
A motivation for such approach can be the difficulty to choose or find an appropriate copula for 
the problem in hand or the fact that the choice of copulas does not inform explicitly about the 
strength of the dependence between the variables involved. Another advantage of the method we 
propose is that we can take tail dependence /extremal coefficients not constant over time and build 
a more realistic time-varying parametric family of models. 

The result shows that the MEV distributions are rich enough to encompass a large variability for 
bivariate dependence one would like to be able to handle. 

2 A multivariate extreme value model 

Let Yi, i = l,...,d > 2. and Zj, j = 1,...,D, be independent and unit Frechet variables and or- , 

d D 

i = l, ..., d, j = 1, ..., D, nonnegative constants. For C constant satisfying C >\J 2_. a \ ■, consider 
the d-dimensional random vector X = (X\, ..., A^) defined as 



d Id 



X t =\/ (afZj) V c - £ <*? \Y t , i = 1, ..., d. 
3=1 \ i=i 



(1) 



We will present in the next result the distribution of the random vector X and its copula 
Cx(wi, ...,Ud) = i 7 x(i r_1 (ui), •■■, F _1 (m,j)), where F^ 1 denotes the generalised inverse of the dis- 
tribution function of A.;, which is a mixture of several complete dependent copulas and total 
independent copulas. 

Proposition 2.1. The mixture model in [w has multivariate extreme value distribution defined 



D ( d 
F x (x u ...,x d ) = l[ /\e" 
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V 3 = 1 J 



(2) 



with its corresponding copula 
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(3) 



Dem. By assuming 1/0 = +oo, we can write from the definition of the model in (fil, for each 
Xi > 0,i = 1, ...,d, 
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F-x(xx,...,Xd) = P 
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(4) 



The assumptions on the distributions of the variables Zj and Yi leads to the expression in d2l 



and, by changing to the variable Ui 



-x~ C 



we obtain the copula function which satisfies the 



max-statbility condition C^_(ui, —,Ud) = Cx(u\, ..., u^), for t > 0. □ 

The tail dependence coefficient between JQ and X^, i,j € {l,...,d} measures the probability 
of occurring extreme values for one random variable given that another assumes an extreme value 
too and is defined as 

Xij = lim P(F Xi (Xi) > u\F X] (Xj) > u), (5) 

where F Xi denotes the distribution function of Xi. In our model, since (Xi,Xj) has bivariate 
extreme value distribution, \j can be computed from the extremal coefficient ejj of (Xi, Xj). This 
is defined by 

Cij = -logF( X . tX .)(x,x) ) x>0, (6) 

where Ft x . x .\ denotes the distribution function of (Xi,Xj). The value of ey doesn't depend on 
x and satisfies e y - = 2 — Xij . 

These dependence measures have extensions for d-dimensional vectors with d > 2 ( |5|10|15|l6] ). 
In the next result we compute Ay for the model in ([l]) and show how such a model for prescribed Ay 
can be constructed. In this way we provide a method for building new multivariate distributions 
and modelling dependence structures. 

Proposition 2.2. (a) The random vector X defined in fn) has tail dependence coefficients 

1 D 

^(^Aaf), s,k€{l,...,d}, s + k. (7) 

i=i 

(b) Given the set of tail dependence coefficients {X s k, s, k € {l,...,d}, s ^ k} of a d-dimensional 
random vector and a constant C > d — 1, there exists a random vector X defined as in Uy with 
proportional tail dependence coefficients {^A s fc, s, k € {1, ..., d}, s ^ k}. 

Dem. First we derive the diagonal of the bivariate copula C s k(u, v) of (X s ,Xk), 1 < s < k < d. 
We have, by taking m = 1 for i £• {s, k} in (J3j, 

D D D D 

(«) A 0) 



A s fc 



C sk (u,u)=u J =1 u J=1 u J=1 =u J=1 (8) 



e fcs = e sfc = - logF (XaA) (C, C) = - logC.^" 1 ^- 1 ) = 2 - i £ (af A oJ. fc) ) , (9) 



Therefore 

which leads to the above presented value for X s k- 

Now, for (b), let A = [X s k]s,k=i,...,d be the symmetric tail dependence coefficients matrix of a 
random vector, where X ss = 1, s = 1, ...,d. Consider a model as in M with a- , i = 1, ...,d, j = 
1, ...,-D = - ~ ' defined as following. Assume that when j <E we don't define the corresponding 
coefficient, n(0) = 0, n(i) — n(i — 1) +d— i, i& {1, ...,d— 1} and n(d) = o~ • ^ e ^ us cons ider 

mi^_j, n ( )_i if j e {ra(0) + 1, ...,n(l)} 



n 



(i) _ I m 2 (5 i _ i: „(i)_2 if j'e{n(l) + l,...,rj(2)} 



j 



(10) 



m,-_i5j_i, n (i_i)_(i_i) if je{rt(i-2) + l,...,n(i-l)}, 



«i- =A iji+ j_ n (i_i), je{n(!-l) + l,..,n(i)}, (11) 



JO 



0, j € i n(i) + 1, 



d(d-l) 



(12) 



with <5i.j = 1 if i = j, Si.j — if i j^ j and m, = y A.; 



i=i+i 



_ d(d-l) 



We remark that n(d — 1) = -j-= — ^ and therefore the set in (13 1 will be empty for i = d — 1. 
Also, since n(0) = 0, the sets in (10 1 will be empty for i = 1. 
For sake of simplicity we can present the coefficients aj in a matrix A = [a, \}i=i,...,d,j=i,...,D as 
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By assuming 5D?__ Oj = if q < p, we can verify that this choice of the coefficients for (jib leads 



D i-l d 

S a i° = Y m i ' + Y Ay < i - 1 + d - i = d - 1 < C, i = 1 

j = l j=l j=i+l 

and, for each pair 1 < s < k < d, 
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+ -^(Xsk A m,) + = — A sfe . 



In fact, as we can seen in the proof, the model can be constructed for any 

d I i—X d 



c>y 5> i+ y, ^ 

i=\ \ j=l j=i+l 



and therefore we have the next two particular cases. 



□ 



Proposition 2.3. (a) Given the set of tail dependence coefficients {X s k, s, k € {1, ...,d}, s ^ k} 

i-l I d \ d 

of a d- dimensional random vector such that >^ \J Xjk + 2_. Ay < 1, for each i = 1, ..., d, 

j=i \fc=j+i / j=i+i 
i/iere exists a random vector X defined as in |7p m£/i ifte same iaiZ dependence coefficients. 

(b) Given the set of tail dependence coefficients {X s k, 8,k € {l,...,d}, s/fc} 0/ a d-dimensional 
random vector such that X s k < 33^ /or eac/i 1 < s < k < d, there exists a random vector X defined 
as in Hy with the same tail dependence coefficients. 



3 Generating MEV distributions 



Example 3.1. Consider the model in UV with a\ 



J 1 ) 



J 1 ) 



2, a? 



4> a 2 



(2) 



2, a[V 



2 = k and (7=1. -By applying the Proposition 1.2, we find the following tail dependence matrix 



A 



1 0.9 0.4 
0.9 1 0.3 
0.4 0.3 1 



for the vector X = (Xi,X2,X3). 
Example 3.2. Lei now 



A 



1 0.2 0.5 0.3 

0.2 1 0.6 0.1 

0.5 0.6 1 0.9 

0.3 0.1 0.9 1 



be the tail dependence coefficients matrix of a A- dimensional random vector. 
We chose 

C = (0.2 + 0.5 + 0.3) V ((0.2 V 0.5 V 0.3) + 0.6 + 0.1) V ((0.2 V 0.5 V 0.3) + (0.6 V 0.1) + 0.9) = 2 

and, as in the proof of part (b) of the Proposition 1.2, we present the coefficients a, , , j = 1, ...,6, 
for Xi along line i of the matrix A, i = 1, ...,4, as follows 



A 



0.2 0.5 0.3 

0.5 0.6 0.1 

0.5 0.6 0.9 

0.5 0.6 0.9 



.,« 



For this choice of a- , i = 1, ...,4, j = 1, ...,6 and C, by applying the Proposition 1.1, we obtain a 
random vector X = {X\, X2, X3, X4) with tail dependence coefficients matrix 



A* = 



1 0.1 0.25 0.15 

0.1 1 0.3 0.05 

0.25 0.3 1 0.45 

0.15 0.05 0.45 1 



Example 3.3. Let now 



A = 



1 0.2 0.1 
0.2 1 0.8 
0.1 0.8 1 



be the tail dependence coefficients matrix of a 3- dimensional random vector. 

If we chose C = (0.2 + 0.1) V ((0.2 V 0.1) + 0.8) = 1, the matrix of coefficients af , i = 1,...,3, j 

1,...,3, 

" 0.2 0.1 
0.2 0.8 
0.2 0.8 

then, by applying Proposition 2.3, the model in U\l still has the given tail dependence coefficients 
matrix A. 



We simulated 1000 observations of X leading to the following bivariate representations. 
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Figure 1: Simulated values for: (Xl,^) (left), (Xl,.^) (center) and (X^^X^) (right). 
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